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1 Introduction 



Exactly solvable problems of quantum mechanics are very interesting and important. Their 
solutions can be found in a straightforward way free of uncertainties caused by perturbation 
methods. The very existence of these solutions is usually connected with symmetries which 
are very attractive subjects by themselves. Moreover, exact solutions present convenient 
bases for expansion of solutions of other problems. 

Another nice property of some quantum mechanical systems is called superintegrability. 
The system is called maximally superintegrable if there exist a sufficient number of alge- 
braically independent operators commuting with the Hamiltonian. This number should be 
equal to In — 2 for the system with n degrees of freedom, and n — 1 of these operators 
should commute amongst each other. 

Many of exactly solvable systems are maximally superintegrable, see, e.g., PQ, and wise 
versa. A perfect example of quantum mechanical system which is both maximally superin- 
tegrable and exactly solvable is the non-relativistic Hydrogen atom (NRHA). The concept 
of superintegrability had been extended to the case of quantum systems which include spin 
[2]. However, the completed relations between the exact solvability and superintegrability 
are not clear yet, especially for the higher order integrals of motion and for the systems 
with n > 2. 

An alternative point of view on the exactly solvable QM systems was created by Gen- 
denshtein [3] who had shown that in many cases such systems admit supersymmetry with 
shape invariance and can be easily solved using tools of supersymmetric QM. In particular 
it is true for the NRHA. 

A consistent maximally superintegrable and exactly solvable problem was discovered by 
Pron'ko and Stroganov [4J. The related quantum mechanical system includes a magnetic 
dipole with spin | (neutron) moving in the field of a stright line current. Like the NRHA, 
this problem admits the hidden symmetry which is more extended then its geometric sym- 
metry pE]. It is supersymmetric also and can be simply solved using the shape invariance 
of its effective potential [5], [§]■ 

The physical and mathematical aspects of the Pron'ko-Stroganov [4J (PS) model was 
multiple investigated. Its experimental realizability was discussed in papers [7] and [SJ, its 
symmetries and supersymmetries was studied in [9] and [10] . enhanced analysis of super- 
symmetric aspects of this model can be found in [llj, etc. Recently the solvable relativistic 
version of the PS model was found [12] which, however, includes a more complicated ex- 
ternal field than in the non-relativistic case. 

Whenever there exist a good model for a particle of spin |, there is a natural desire to 
generalize it for higher spins. In the case of the PS model the time needed to realize this 
idea was thirty long years. The main hidden dangers in this generalization were connected 
with the fact that the direct change of the Pauli interaction term present in the PS model 
by the analogous term for higher spin leads to the system which is not exactly solvable and 
its finite trajectories are not closed [13]. 

In paper [14] exactly solvable generalizations of the PS model to the case of arbitrary 
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spin have been formulated. The price paid for this progress was the essential complication 
of the Pauli interaction term present in the PS model. However, there are good physical 
arguments for such a complication in the case of spin s > ~, in both the non-relativistic 
[H] and relativistic [15] approaches. 

A natural question arises whether the models proposed in [T4] can be effectively inte- 
grated using the tools of SUSY quantum mechanics, like it has been done [6] for the PS 
model for neutral particle of spin |. This question is especially provocative for us since in 
the recent papers [16] and [IT] an effective classification of shape invariant matrix potentials 
was presented. Thus if the matrix potentials found in [H] are shape invariant they should 
be nothing but particular cases obtained by this classification. 

In the present paper the supersymmetric aspects of the maximally superintegrable sys- 
tems for arbitrary spin (proposed in [14]) are discussed. It is shown that these systems 
are shape invariant and can be easily integrated using the tools of SUSY QM. Using the 
results of paper [16] we construct exact solutions of the mentioned systems directly for ar- 
bitrary spin. In addition, a more straightforward and refined formulation of these systems 
is proposed together with the proper physical interpretation. 



2 Hidden symmetry of the PS model 

The PS model is based on the following version of the Schrodinger-Pauli Hamiltonian: 



H = Pl±A + X SxV o SyX (2.1) 
2m H 



which is the Hamiltonian of a neutral spinor anomalously interacting with the magnetic 
field generated by a straight line current directed along the z coordinate axis. Here p x = 
—i-^,Py = — ij^, r 2 = x 2 + y 2 ] S x = \o\ and S y = \ai are matrices of spin |, a± and oi 
are Pauli matrices, A is the integrated coupling constant. 

The last term in ( 12. ip is the Pauli interaction term AS • H where the magnetic field H 
has the following components which we write ignoring the constant multiplier included into 
the parameter A: 

H x ~ H y ~ H z = 0. (2.2) 

Hamiltonian ( 12. ip is invariant w.r.t. rotations around the z-axis since it commutes with 
the z component of the total angular momentum 

J z = X1P2 - X2P1 + S z (2.3) 
where S z = -03. In addition, it admits two more constants of motion [4j, namely: 

I Tfl I Tfl 

A x = -(JzPx + PxJz) + —n(n)y, A y = -(J z p y + p y J z ) - — /x(n)x (2.4) 
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where ;u(n) = X(S x n y — S y n x ), n x = -,n y = -, and n = (n x , n y ). Operators ( 12. 3p and (12 ,4p 
commute with % and satisfy the following commutation relations: 

\Jzi A x ] = iA y , [J z , A y ] = —iA x , [A x ,A y ] = —J Z T-L. (2-5) 

In other words, operators (12.41) and (12.51) represent algebra o(l,2) on the spaces of eigen- 
vectors of H. and cause the degeneration of the Hamiltonian eigenvalues. And this hidden 
symmetry with respect to algebra o(l,2) makes the PS model maximally superintegrable 
and exactly solvable. 

It is interesting to note that up to unitary equivalence operator (13. ip is the only plane 
Hamiltonian for neutral particle of spin | interacting with an external field, which commutes 
with J3 and admits a hidden symmetry w.r.t. algebra o(l,2). This statement in fact was 
proven in paper [4J where, however, a more general Hamiltonian has been proposed, namely: 

-,2 1 j 



u = Px+Py + a S x yS y x + b S x x + S y y 
2m r 1 r l 



Operator (I2.6P does commute with (12 .4p and (12. 5p and includes two arbitrary parameters 
a and b. But it is unitary equivalent to (13. ip since 

U = UUU^ with U = cosfl + LS.sinfl (2.7) 

where S z = and 9 is a real parameter such that a = X cos 9 and b = Asin#. 

3 SUSY aspects of PS model 

Let us apply transformation (12 ,7p . For our purposes it is convenient to set 9 = tt/2, then 

n = pJ±A + x ^±M. ( 3.i) 

2m r 2 

Then we consider the eigenvalue problem for Hamiltonian ( 13. ip : 

Ui\) = Sip (3.2) 

where ip = ip(x,y) is a two-component wave function. This function is supposed to be 
square integrable and vanish at x = y = 0. 
Introducing the polar coordinates 

x = rcos9, y = rsin9 (3-3) 

and expanding ip via eigenf unctions of the angular momentum operator J z : 

,1, -cui, ii, - — ( e Mi("- 3W1 ^ (oa) 
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where C k are constants, 0i and 02 are functions of r and summation is imposed over the 
repeated indices k = ±|,±|,±|---, equation ( 13 .2p is reduced to the following system of 
decoupled equations for ip K : 

(Q2 1 X\ ~ 

where £ = 2mS, A = 2mA. 

Equation (13. ip and its reduced version (I3.5P are invariant with respect to the space 
reflection x — > x,y — > —y provided ip{ x ,y) an d ipk{r) cotransform as 

ip(x, y) Aip(x, -y), ip k (r) Atp_ k (r) (3.6) 

where A = <j\. 

Transformations (I3.6P commute with Hamiltonians % and T-L K , anticommute with J 3 
and are invertible. Thus we can restrict ourselves to equation (13.51) for positive half integer 
k. Then solutions for negative k can be obtained using (13.61) . 

Hamiltonian T-t K can be factorized as 

U K = a+a K + c K (3.7) 

where 

«. = ^ + W„ < = ~ + W„ C = -p^ (3-8) 
and W is the matrix superpotential 

W K = -U 3 - - -fc^. (3.9) 

2r 2k + 1 r 

One more important property of 1-L R is its shape invariance, i.e., its superpartner is 
equal to T-L K+ \ up to a constant term: 

d 2 1 A 

Ht = a-a+ + c K = -— + + 1)(k + 1 - a 3 )- + a x - = V, K+l + C K 

where C K = c K — c K+ x. Thus equation (13.51) can be solved using the standard tools of SSQM. 
Namely, the ground state vector is defined as a square integrable solution of the first order 
equation 

a-iP Kt0 {r) =(^ + WA ^,o(r) = (3.10) 

and thus it has the following components 

fa = u^K^u), 2 = u K+1 K (u) (3.11) 
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Xr 



X 2 



(2fc+l) 2 ' 



where Kn and K\ are the modified Bessel functions, u — , 

In view of (13 .7p function V«,o(0 solves also equation (13. 5 j) with £ = c K 
Solutions which correspond to n th exited state can be represented as 

AA r ) = at a t+i ■ ■ ■ <4+n-iVw,o(z)- (3.12) 
The corresponding eigenvalue S n is given by the following formula: 

n-l r 2 

£ n = = - (2K + 2n + 1)2 - (3-13) 



i=0 



4 Integrable models for neutral vector bosons 

The PS model considered in the previous section describes the interaction of a neutral 
particle of spin | with the field of straight line current. How we can generalize Hamiltonian 
(13. ip (or the initial Hamiltonian (12. ip which is unitary equivalent to (13. ip ) to the case of 
higher spins? The standard idea is simple to change in (13.11) matrices S a = \a a by matrices 
of higher spin. However, in this way we obtain a model which is not either exactly solvable 
[13J or maximally superintegrable. 

To obtain a maximally superintegrable analogues of (12. ip for higher spins it is necessary 
to make a rather non-trivial generalization of the Pauli interaction term in the Hamiltonian 
(13. ip . including multipole magnetic interactions [14] . 

In this section we present an analogue of the PS model for particle of spin 1. This 
model is both maximally super integrable and shape invariant. It is based on the following 
Hamiltonian 

K = ^U-^(n) (4-1) 
2m r 

where 

fi s (n) = ^(n) = /i(2(S x n) 2 - 1) + A(2(S • n) 2 - 1). (4.2) 

Here fi and A are arbitrary real parameters, S • n = Sini + 62^2 and S x n = Sin 2 — S^ni, 
5*i and S 2 are matrices of spin 1: 

S y = -j=\l -1 I , S z =\ I . (4.3) 

Just Hamiltonian (14. ip . (14. 2p is the absolutely super integrable analogue of (12. ip for 
spin one. Matrix ( I4.2p depends on n = (-, ^) and satisfies the following conditions: 

Mn),J 2 ] = 0, (4.4) 
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Hs(n)S z + S z /i s (n) = 0. 



(4.5) 



Thus Hamiltonian (14. ip commutes with operators (12. 3p and (12. 4p where /i(n) — >■ /xi(n) and 
S z is the matrix given in ( 14. 3p . So this Hamiltonian admits dynamical symmetry w.r.t. 
algebra o(l,2). 

It is possible to show that, up to unitary equivalence, equations (14. ip . (14. 2p represent 
the most general form of plane Hamiltonian for neutral particle of spin 1, which admits this 
symmetry, see section 7. We will see that, in addition, this Hamiltonian is shape invariant. 
Let us discuss its physical content. 

The physical sense of the PS model for spin | is absolutely clear. The related Hamilto- 
nian (I3.2p includes the Pauli term ~ S • H and corresponds to the interaction of a neutral 
particle (having a non-trivial dipole moment) with the field of the constant and straight 
line current. 

Physical content of superintegrable models for higher spins is much more sophisticated. 
The interaction term of Hamiltonian (I3.2p has nothing to do with the Pauli term and needs 
another interpretation. 

Let us consider in more detail the case /i = 0, A = u > 0. The corresponding Hamilto- 
nian (14. ip takes the following form: 

Hi = A^ +tt (Wj£_l). (4 . 6) 



2m \ r 3 r t 

and can be represented as 

n, V 2 Xl +vl 2 . „ dE a u . 
n ' = —2^ + uQ ^ + 3 dlvE (4 - 7) 

where 

2 

Qab = S a S b + S b S a - -s(s + l)S ab 
is the tensor of quadruple interaction, 

E=(^,0), (4.8) 

and the temporary notations x — Xx, y — Xi are used. 

In accordance with (14.71) Tii can be interpreted as a Hamiltonian of spin-one particle 
which has neither minimal nor dipole interaction with the external field. However, this 
particle is supposed to have a quadruple and Darwin interaction, the latest is represented 
by the last term in (14.71) . 

A more delicate question is related to the physical realizability of the vector field (14.81) 
included into the interaction terms. If we suppose that it is the classical Maxwell elec- 
tric field, the corresponding charge density should be proportional to divE = -. Such 
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charge density is seemed to be hardly realized experimentally. However, vector (14.81) per- 
fectly solves field equations of generalized Maxwell electrodynamics modified by presence 
of Chern-Simons term, and also equations of axion electrodynamics with trivial current 
and charge densities [T8] . 

Let us present one more expression of the interaction term via physical fields. Setting 
A = 0, fi = oj and restricting ourselves to /i > we can rewrite Hamiltonian (14.1 1) in the 
following form: 



4 _Pl+P 2 y , ,2(S-H) 2 -H 2 

ft i — — h A- 



2m 



H 



(4.9) 



where H is the vector of magnetic field whose components are defined in equation (I2.2p . 
Notice that Hamiltonians (14. 6p and ( I4.9P are unitary equivalent. Namely, 



(4.10) 



where s — 1. 

The last term in (14.91) represents a nonlinear interaction with the well defined external 
magnetic field (12. 2p . This field solves the Maxwell equations with a constant straight line 
current. 

Let us remind that a nonlinear generalization of the Pauli interaction is requested also 
in the relativistic description of spin-one particle interacting with the constant magnetic 
field 



5 Absolutely super-integrable models for spin | 



Let us discuss an exactly solvable model for particle of spin |. The corresponding Hamil- 



tonian has the generic form (14. ip with 

^(n) = (u + iiSl) (7S x n - 4(S x n) 3 ) . 

Here /i and v are arbitrary parameters, S x , S y and S z are the 4x4 matrices of spin 
can be chosen in the following form: 



(5.1) 
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Matrix (15.1 j) satisfies conditions (14. 4p and (14. 5p which are necessary and sufficient for 
existence of constants of motion (I2.3P and (12 .4p [13]. Thus the system whose Hamiltonian 
is given by equations (14. ip and (15. ip actually is maximally superintegrable. 

The interaction term of the considered Hamiltonian can be represented in terms of 
external fields. In the case = 0, v = | ^ we have: 

where to = H is the vector of magnetic field whose components are defined in equation 
(12. 2p and S is the spin | vector with components (15.21) . 

In addition to the standard Pauli term uS ■ H Hamiltonian (I5.3P includes the additional 

(S H)^ 

interaction ~ jjr - which is non-linear in magnetic field. 

Alternatively, for /i = — |z/ we obtain the following representation: 



where 




is the octuple interaction tensor (the summation is imposed over all possible permutations 
of indices a, b and c), and H = (x 2 lnr, — X\ lnr, 0) . 

The vector field H perfectly solves equations of the axion electrodynamics [18] . On the 
other hand, being treated as the classical magnetic field it requires a current whose density 
grows algorithmically with growing of r. 



6 Superintegrable systems for arbitrary spin 

In paper [2] the generalization of Hamiltonian (12. ip for arbitrary spin s was proposed, 
which admit dynamical symmetry w.r.t. algebra o(l,2). This Hamiltonian has the generic 
form (14. ip where fi s (n) is a (2s + 1) x (2s + 1) dimensional matrix depending on n = (-, ^) 
and satisfying conditions (I4.4p . (14. 5p with S z being the z component of spin s vector, i.e., 
the matrix 

S z = diag(s, s — 1, s — 2, ... — s). (6.1) 

It can be verified by direct calculations that if conditions (14. 4p and (14 .5p are satisfied 
then Hamiltonian (14. ip commutes with operators (I2.3p and (12.41) (where S z is matrix (16. ip 
and yu(n) — > /i s (n)), and these operators do satisfy relations (12.51) . 

We shall refine and extend the results of paper [13]. First, we present a straightforward 
formulation of Hamiltonians (14. ip for arbitrary spin. Secondly, we prove that the number of 
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arbitrary parameters present in the models found in [2] can be effectively reduced applying 
a unitary transformation. At the third place the shape invariance of Hamiltonians (14. ip 
will be proven and the corresponding eigenvalue problems will be solved algebraically. 

Let us find Hamiltonians (14. ip for arbitrary spin. The condition ( 14. 4 j) is satisfied iff 
/i s (n) is a function of S • n = S x n x + S y n y , S x n = S x n y — S y n x and S z . At the first step 
we restrict ourselves to matrices /i s (n) which are polynomials in S • n. 

It is convenient to represent matrix // s (n) in the following form: 



[n) = C » A » (6-2) 



where c u are unknown coefficients and A u are projectors onto the eigenspaces of matrix 
S ■ n corresponding to the eigenvalue v (v, v' — s, s — 1, ... — s): 

* IP" (6.3) 

Matrix (JO} by construction satisfies condition ( I4.4p . Substituting (16. 2 p into ( 14. 5 p and 
using the identities 

S z fi s (n) + /i a (n)S z = 2S Z /J, s (n) + [^(n), S z ) 

and 



[A v , S z ] = -S Z {2A U - A v+1 - A^i) + ^ (n x S y - n y S x ) (A u+1 - A v _i) (6.4) 

we easily find the following condition for coefficients c u (note that projectors A u are orthog- 
onal and linearly independent): 

Cv = -Cv-i, v = s,s - 1, 1 - s. (6.5) 

In accordance with (16. 2p and (16.51) the general expression for /i s (n) can be given by the 
following equation: 

^ s (n) = A^(-l)MA, (6.6) 

V 

where [u] is the entire part of i>, v = s, s — 1, — s. 

Equations (14. ip . (16. 6p present Hamiltonians for arbitrary spin s admitting dynamical 
symmetry w.r.t. algebra o(l,2). Potentials ( 16. 6 p are defined up to arbitrary parameter A 
which can be associated with the coupling constant. 



2 Formula (|6.4p is a particular case of equation (A. 2) in [TO] corresponding to ^ = n x , ^ = n y , ^ = 
0, 5*23 = S x ,Ss\ = S y . This formula is also a consequence of equation (7.18) from book [3D] 
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However in accordance with the results of paper [J3] the general solution of the deter- 
mining equations (I4.4p . (14. 5p depends on 2s + 1 arbitrary real parameters. This statement 
is in accordance with the fact that there exist exactly 2s + 1 linearly independent matrices 
anticommuting with S z . 

But our analysis admits a straightforward and simple extension to the generic case. 
Indeed, equations (I4.4p . (14. 5p are invariant w.r.t. multiplying /i s (n) by arbitrary power of 
matrix S z . Thus, starting with (16. 6p we can immediately construct 2s + 1-parametrical 
solutions of this equation: 

/} s (n) = ^(b v B v + id v C v )ii s (n) (6.7) 

where /i s (n) are matrices (16. 6p . b u and d u are arbitrary (real) parameters, B v and C u are 
projectors, polynomial in S z : 

Q2 _ 12 o o2 _ 12 

^=Il^-^f> ^ = -11^—^. = 8,8-1,8-2,..., U,U'>0. 

- LJ - u z — v v - LJ - v A — V A 

Notice that C u includes only odd powers of S z and so they anticommute with yU s (n). Thus 
matrix /2 s (n) is hermitian. 

Using the orthogonality relations 

it is easy to prove that all terms in sums (16. 7p are linearly independent. Thus matrix fi s (n) 
includes exactly 2s + 1 essential parameters and so it should be equivalent to the analogous 
matrix found in [13] , see equation (29) there. However the number of these parameters can 
be reduced using the unitary transformation 

fl s (n) ->• UpL s {rv)U ] = KB v fi a (p) = P>,(p) (6-8) 

u>0 



where \ u = y^bl + d 2 v and U = ^ v>Q (cos9 v B v + i$m6 u C v ). Taking into account that // s (n) 
commutes with B v and anticommutes with C v it is easy to make sure that relation (16.81) is 
true provided parameters 9 V satisfy the conditions b u = X v cos 9 V and d u = X u sm9 u . 

It follows from (16. 8p that, up to unitary equivalence, matrix /t s (n) includes only s + | 
or s + 1 arbitrary parameters for half integer or integer spin respectively. A particular case 
of this statement has been proven in section 2, see equations (12. 6p and (12. 7p . 

Thus the Hamiltonian for superintegrable model of arbitrary spin can be represented in 
the following form: 

H s = 4^ + -^(n) (6.9) 
2m r 
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where /2 s (n) is the matrix defined by equations (16.81) and (16. 6p . Hamiltonian (16 .9p is invari- 
ant w.r.t. algebra o(2,l) whose basis elements are given by equations (I2.3P and (12 .4p with 
/x(n) fi s (n). 

Notice that in particular case s = 1 Hamiltonian (16. 9p is reduced to the operator defined 
by equations (14.11) and (14.21) where /i = |(Ao — Ai) and A = |(Ao + Ai). For s = 3/2 matrix 
(16. 8 p can be transformed to the form (15. ip with ji = |(A| — Ai ), v = |Ai — ^A| by applying 
transformation (14.101) . 



7 Dual shape invariance and exact solutions 

Let us consider the eigenvalue problem (13. 2 p for Hamiltonian (14. ip . (16. 6p . Like in section 2 
we again introduce polar coordinates (13. 3 p and expand solutions via eigenvectors of operator 
J z (I2.3p where S z is the matrix of spin s defined by equation (16.11) . These eigenvectors can 
be represented as: 

^ = 4=exp(i(«-S;)0)$«(r) (7.1) 



where 

$ re (r) = coulomb(0 s , S _ X , (7.2) 

Substituting (16. 9p . (17. ip . (17. 2p into (I3.2p we obtain the following equation for radial 
functions 

= (- + V^j $ K = £<b K (7.3) 
where 8 = and 

V K = i(k - S z ) 2 - - ) — + /i s -, /i s = 2m/i s (n)| n - . (7.4) 
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Both the initial equation (13. ip and equation (17. 3p are invariant w.r.t. the space reflection 
transformation (13.61) where A = /i s (n)/2mA. Thus it is reasonable to search for solutions of 
(17. 3p for non-negative k since solutions for negative k can be obtained using transformation 
(13751). 

Hamiltonian H K and matrix S 2 commute between themselves and so they have a mutual 
system of eigenfunctions. Matrix S 2 is diagonal and its eigenf unctions ifi u corresponding to 
eigenvalues v 2 = s 2 , (s — l) 2 , (s — 2) 2 , ... have two or one non-zero components for v 2 > 
and v = respectively. In the standard representation of spin matrix with diagonal S z 
and symmetric , matrix fl s is symmetric and antidiagonal, and has the unit non-zero 
entries. 



3 See, e.g., QH, eq. (22), or [3D] eq.(4.65) 
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Thus the eigenvalue problem (17. 3p can be decoupled to the following equations: 
/ d 2 \ 

Uk^w = —z-z + Vk 4>k,v = £^k, v (7.5) 



dr 

where ij) KV = I J v ] are two-component functions, and 



(k — w 3 ) 2 — 7 A , 
V«r = - r ~ + -^ ( 7 - 6 ) 



For v = the corresponding reduced potential is one-dimensional: 
^ ) 2 ~ 1 + A 



Potentials (17. 6p were discussed in paper [16] where parameter i/ was denoted as fi + \. 
These potentials are shape invariant and can be expressed via superpotentials as follows: 

Vk,u = Wl tV - W' KtV + c K (7.8) 

where 

and c K is the parameter given in equation (13. 8p . 

Shape invariance of can be easily proven since, in addition to the representation 
(I7.8P the following equation holds true: 

K, = Wl„ + W'^ = V k+h „ + C K (7.10) 

where C R = c k+1 - c K . 

One dimensional potential (17. Tp is shape invariant too. It can be represented in the 
standard form (17.81) with 

It is important to note that potential ( 17. 6p is invariant w.r.t. the change v — y K, K — > i>, 
while superpotential (17. 9p is not invariant w.r.t. this change, since 

^=^3- * , (7.12) 
r 2v + 1 2r 

Thus, in addition to (17. 8p . there exist the alternative representations of potential (17. 6p 
via superpotential, i.e., equation ( I7.8P where W KtU is changed by superpotential (17. 12ft : 

V^ = Wl K -W'^ K + c ul = - * (7.13) 

12 



In other words, potential ( l7.6|) appears to be shape invariant w.r.t. the shifts of two 
parameters, i.e., k and v. This is a particular case of the dual shape invariance phenomena 
discovered in [T6] . 

Using representations (I7.8P and (I7.13P we easy find the ground state vectors of Hamil- 
tonian (17. 5p . Solving equations (I3.10p with superpotentials (17. 9p and (I7.12p we obtain the 
following components of the ground state vectors i\P v = column(0°, 0°_ v ), v = s, s — 1, s — 



2,...,u> 0: 



Ar \ . , i +1 / Ar 



(7.14) 



for superpotential (17.91) . and 

2 \ 2u+ 1 r 




(7.15) 



for superpotential (I7.12p . were d u are integration constants. Solution for v = 0, i.e., the 
component 0°, is given by the following equation: 

0[J = du K+ ^ exp (-Ar) , k = 0,1,2,... (7.16) 



Functions (I7.14p and (I7.15P are square integrable for A > and arbitrary integer or 
half integer k > 0. The same is true for (I7.16p . However, functions (I7.14p for k < v and 
functions (17.151) for k > v do not vanish at r = 0. So such values of parameters v and k 
should be excluded, as it is indicated in the r.h.s. of the discussed equations. 

Vectors for exited states and the corresponding energy levels again are given by relations 
(I3.12p . (13. 8p and (13 . 13[) respectively where W K should be replaced by superpotential (17. 9p 
or (I7.12p . In accordance with the analysis presented in [16] all such vectors are square 
integrable and vanish at r = 0. 

Thus we find the eigenvectors of Hamiltonian (14. ip . (I6.6P which are given by equations 
(17. ip . (17. 2p and (I7.14p -( f?.16p . The corresponding solutions for negative k can be easily 
found using transformation (13.61) where A = /i s (n). 

In complete analogy with the above we can solve the eigenvalue problem for the more 
general Hamiltonian (16. 9p including [s + 1] arbitrary parameters \ u . Actually, to this effect 
it is sufficient to change A — > \ u in all formulae (17. 6p . (I7.9p -( I7TT5]) and (I3.13p . 

8 Discussion 

Shape invariant matrix potentials classified in papers [16] and [17] can appear in many 
realistic integrable models of quantum mechanics. In the present paper we apply the 
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results of these papers to the interesting class of exactly solvable systems found in [T4] . 
These systems are maximally superintegrable and generalize the PS model to the case of 
arbitrary spin. 

We find it be interesting to study these long-awaited generalized models in more detail. 
In particular, to search for their solutions, to examine their consistency, and to verify 
weather the supersymmetry and shape invariance of the spin ~ model are kept in the case 
of arbitrary spin. In addition, it is important to understand the physical content of the 
models proposed in [T4] . Just the tasks enumerated above are the subjects of the present 
paper. 

First we refine the results of paper [TJ] taking into account equivalence relations w.r.t. 
unitary transformations realized by constant matrices. It is shown that up to such equiv- 
alence the superintegrable models for arbitrary spin s include s + 1 or s + | arbitrary 
parameters for integer or half integer spins correspondingly, while in [T4] these models 
include 2s + 1 parameter. 

The cases s = 1 and s = | are considered in more detail. We present the corre- 
sponding Hamiltonians in the form which is convenient for physical interpretation. This 
interpretation is proposed using two alternative ways. First, it is possible to treat f l4.ll) as 
a Hamiltonian of spin-one particle with zero charge and zero dipole momentum, which has 
a non-trivial quadruple interaction with the external field. Another possibility is to repre- 
sent the Hamiltonian in form (I4.9p including the interaction term nonlinear in the external 
magnetic field. Analogous alternatives exist for the models for spin | particle considered 
in section |5]and also for arbitrary spin, i.e., it is possible to interpret the related potentials 
as results of either multipole or non-linear interaction of a neutral particle with an external 
field. 

It is shown that the considered models are shape invariant. Moreover, their effective 
potentials appear to be particular cases of the matrix shape invariant potentials classified 
in [IB]. Using this fact it is possible to construct exact solutions of these models in a 
simple and straightforward way. This program has been realized in section [7] immediately 
for arbitrary spin. 

A specific property of potentials (17.61) is their dual shape invariance, i.e., the existence 
of two non-equivalent corresponding superpotentials. Exactly this property enables to 
find good solutions for all combinations of eigenvalues of the spin and orbital momentum 
operators. 

Notice that supersymmetric and superconformal aspects of various planar systems with 
arbitrary spin were intensively discussed in literature, see, e.g., |2T] and [22]. We consider 
special class of such systems which are chargeless, have non-trivial multipole moments, and 
can be integrated in closed form using their supersymmetry and shape invariance. 

In the present paper the shape invariant matrix potentials classified in [IB] and [FT] 
are used to solve explicitly the countable set of maximally superintegrable models. In fact 
the number of exactly solvable models with matrix superpotentials some of which have 
the clear physical significance is much more extended. Among them there are also the 
relativistic models discussed in [12] and many others. We plane to present such models in 
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the following publications. 
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